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ABSTRACT 


An  approximate  analytical  solution  is  obtained  for  hypersonic  flow 
past  a slender  elliptic  cone  at  small  angle  of  attack.  The  analysis  is 
based  on  perturbations  of  hypersonic  flow  past  a circular  cone  aligned 
with  the  free  stream,  the  perturbations  stemming  from  the  linear  combina- 
tion of  small  angle  of  attack  and  small  cross-section  eccentricity.  By 
means  of  previously  obtained  hypersonic  approximations  for  the  basic- 
cone  problem,  closed-form  approximate  solutions  for  the  perturbation 
equations  are  obtained  within  the  framework  of  hypersonic  small -disturbance 
theory.  Results  for  the  shock  shape,  shock-layer  structure,  and  surface 
conditions  are  presented,  together  with  comparisons  with  experimental 
data. 


□ □ 


ACKNOWLEDGEMENTS 


This  investigation  was  sponsored  by  the  United  States  Air  Force, 
Office  of  Scientific  Research,  under  Grant  No.  AFOSR  77-3468.  This  report 
is  a slightly  modified  version  of  the  Master  of  Science  thesis  by  Mr. 
Hsiung  Ming  Lee.  The  authors  thank  Professor  Martin  Jischke  for  numerous 
valuable  discussions.  Thanks  go  also  to  our  project  monitor.  Dr.  Donald 
Daniel,  for  his  encouragement  and  comments. 


TABLE  OF  CONTENTS 


ACKNOWLEDGEMENTS  iii 

LIST  OF  ILLUSTRATIONS vi 

LIST  OF  SYMBOLS vi  i 

1.  INTRODUCTION  1 

2.  UNYAWED  ELLIPTIC  CONE  BODY 

AND  SHOCK  GEOMETRY  4 

2.1  Body  Geometry 4 

2.2  Shock  Geometry  7 

3.  BOUNDARY  CONDITIONS  9 

3.1  Expansions  for  the  Flow  Variables 9 

3.2  Free  Stream  Normal  Velocity 

at  the  Shock  9 

3.3  Downstream  Pressure  at  the  Shock 10 

3.4  Downstream  Density  at  the  Shock  11 

3.5  Downstream  Velocity  at  the  Shock 11 

4.  PERTURBATION  EQUATIONS  14 

4.1  Energy  Equation 14 

4.2  Momentum  Equation  16 

4.3  Continuity  Equation  18 

5.  APPROXIMATE  SOLUTION  FOR  HYPERSONIC 

FLOW 23 

5.1  Basic  Cone-Flow  Approximations 23 

5.2  Evaluations  of  Integrals 23 

5.3  Approximation  for  the  Velocity 24 

5.4  Shock  Eccentricity  Factor  26 

5.5  Shock-Layer  Velocity  Profiles  27 

5.6  Evaluation  of  the  Approximate  Analysis 33 

6.  SURFACE  PRESSURE  ON  THE  ELLIPTIC  CONE 38 

6.1  Surface  Perturbation  Pressure 

Coefficient 38 

6.2  Comparison  with  Experiment 41 

6.3  Drag  on  the  Elliptic  Cone 43 

7.  ELLIPTIC  CONE  AT  ANGLE  OF  ATTACK 47 

7.1  Superposition  of  Results 47 

' l . 7.2  Pressure  on  the  Body  Surface 48 

7.3  Comparison  with  Other  Results  49 

7.4  Normal  Force  on  the  Elliptic  Cone 54 


iv 


j 


Figure 


LIST  OF  ILLUSTRATIONS 


Cone  Coordinates  and  Geometry 

Two-Term  Fourier  Representation  of  the  Ellipse 

Elliptic-Cone  Shock  Eccentricity  Factor, 
y 3 1 .4 

Radial  Perturbation  Velocity,  y = 1.4 

Polar  Perturbation  Velocity  Components,  y=  1.4 

Azimuthal  Perturbation  Velocity  Component,  y * 1.4 

Azimuthal  Velocity  at  the  Body  Surface, y=  1-4 

Correction  Velocity  Ratio,  y = 1.4 

Perturbation  Pressure  Coefficient  on  the  Body  Surface 

Pressure  Coefficient  on  Body  Surface, 

Model  I,  « 6,  y 3 1.4 

Pressure  Coefficient  or.  Body  Surface, 

Model  II,  M = 6,  y = 1.4 

oo 

Pressure  Coefficient  on  Body  Surface, 

Model  I,  M = 3.09,  y - 1.4 
00 

Pressure  Coefficient  on  Body  Surface, 

Model  II,  M = 3.09,  y = 1.4 
00 

Angle-of -Attack  Perturbation  Surface  Pressure  Coefficient 
Y = 1 .4 

Comparison  of  Surface  Pressure  with  Martellucci's 
Results,  Model  I,  M = 6,  y = 1*4 

oo 

Comparison  of  Surface  Pressure  with  Martellucci's 
Results,  Model  II, M =6,y=1-4 


r- 


< 


LIST  OF  SYMBOLS 

a tangent  of  semivertex  angle  of  semiminor  axis 

&o  zeroth  order  isentropic  speed  of  sound 

A basic  area  of  finite  cone 

b tangent  of  semivertex  angle  of  semimajor  axis 

moment  coefficient 
normal  force  coefficient 
cy  constant  density  specific  neat 

c surface  pressure  coefficient 

Ar 

Cp  surface  pressure  coefficient  associated  with  angle  of  attack 

D drag  force 

e eccentricity  of  elliptic  cone 

F force  on  finite  cone 

F,  constant,  see  Equation  (4.10) 

'X/ 

F.j  constant,  see  Equation  (7.6) 

constant,  see  Equation  (5.10) 

g shock  eccentricity  factor  associated  with  elliptic  cone 

g shock  eccentricity  factor  associated  with  angle  of  attack 

Hq( 9 ) function,  see  Equation  (4.32a) 

H^(e)  function,  see  Equation  (4.32b) 

r ( e ) function  of  integration,  see  Equation  (4.22) 

Kj  hypersonic  similarity  parameter 

M Mach  number 

M moment  about  the  cone  vertex 

N constant,  see  Equation  (5.24) 


-0 


vi  i 


n 


unit  normal  vector 


' 


p pressure 

p pressure  associated  with  angle  of  attack 

r spherical  coordinate 

* 

s entropy 

T temperature 

u,v,w  spherical  velocity  components 

u,v,w  spherical  velocity  components  associated  with  angle  of  attack 
V velocity 

X(e)  function,  see  Equation  (4.33) 
x,y,z  Cartesian  coordinate 

Greek  Letters 

a angle  of  attack 

8 unperturbed  shock  angle 

y ratio  of  specific  heats 

5 semivertex  angle  of  basic  circular  cone 

e measure  of  eccentricity 

0,4)  spherical  coordinates 

5 p^/pU)  > density  ratio  across  the  shock 
p density 

a 6/5 

n vorticity 


viii 


APPROXIMATION  FOR  HYPERSONIC  FLOW 
PAST  AN  ELLIPTIC  CONE 


1 . INTRODUCTION 

For  supersonic  flows  past  bodies  without  axial  symmetry,  the  ellip- 
tic cone  is  a basic  body  shape.  Its  counterpart,  the  circular  cone,  is 
a basic  axisymmetric  body  shape  with  flow  field  properties  that  are  exten- 
sively tabulated.  Even  the  supersonic  flow  past  a circular  cone  at  small 
angle  of  attack  is  extensively  tabulated  and  fairly  well  understood.  On 
the  other  hand,  the  properties  of  the  supersonic  flow  field  past  an 
elliptic  cone  are  not  extensively  tabulated,  at  least  in  comparison  with 
the  circular  cone  at  angle  of  attack.  Although  numerous  papers  have  been 
directed  towards  supersonic  flows  past  elliptic  cones,  their  goals  have 
been  specific,  and  no  general  or  comprehensive  flow  field  calculations 
have  been  set  forth.  The  ourpose  of  this  work  is  to  partially  remedy 
this  situation  and  to  present  an  approximate  analytic  solution  that 
illustrates  the  general  flow  field  features  of  hypersonic  flow  past  a 
slender  elliptic  cone  with  small  eccentricity. 

Supersonic  flow  past  slender  cones  with  arbitrary  cross  section  can 
be  treated  with  some  generality  by  means  of  linearized  theory  [1],  and 
linearized  theory  can  be  extended  to  second  order  [2]  to  account  some- 
what for  weak  nonlinear  effects.  These  results,  however,  are  frequently 
not  appropriate  for  the  hypersonic  Mach  numbers  and  flow  deflections  of 
practical  concern. 

Another  approach  to  the  problem  deals  with  flows  that  deviate  slightly 
from  the  basic  axisymmetric  flow  past  a circular  cone.  The  so-called 


-2- 


scheme  of  linearized  characteristics  was  applied  by  Ferri , Ness,  and 
Kaplita  [3]  to  several  conical  bodies  with  non-axi symmetric  cross  sec- 
tions. This  method  is  subject  to  a number  of  criticisms  for  cross-sec- 
tion areas  that  deviate  significantly  from  circles  [4],  and  thus  a modi- 
fication of  the  linearized  characteristics  method  was  applied  to  elliptic 
cones  at  angle  of  attack  by  Martel lucci  [5],  In  the  above  use  of  the 
linearized  characteristics  method,  the  perturbation  equations  were 
solved  numerically.  Chapkis  [6],  applying  the  linearized  characteristics 
method,  used  hypersonic  approximations  for  the  basic  cone  flow  to  obtain 
relatively  simple  specific  results  for  an  elliptic  cone. 

Besides  the  above  methods  of  computation,  there  are  numerical  schemes 
for  integrating  the  complete  governing  gas  dynamic  equations.  Two  notable 
schemes  applied  to  elliptic  cones  are  those  of  Stocker  and  Mauger  [7] 
and  Babenko  et  al . [8] . There  are  also  semi -empirical  methods  for  deal- 
ing with  certain  features  of  supersonic  flows  past  elliptic  cones,  such 


as  tangent-cone  methods,  equivalent  circular-cone  methods,  and  the  method 
of  Kaattari  [9].  Whichever  method  has  been  utilized  to  date,  the  general 
features  of  the  supersonic  flow  field  past  an  elliptic  cone,  showing 
effects  of  Mach  number,  cone  angle,  and  ellipse  eccentricity  on  the  shock 
shape,  the  shock-layer  structure,  and  the  surface  condi tions,  have  not 
been  delineated. 

In  this  undertaking,  we  wish  to  start  with  the  small -perturbation 
equations  for  perturbed  flow  past  a basic  circular  cone  at  zero  angle 
of  attack  and  take  the  approach  of  Chapkis  [6]  for  the  elliptic  cone. 

We  shall  use  improved  approximations  for  the  basic  cone  flow  [10,11]  and 
obtain  approximate  analytical  solutions  for  the  perturbation  equations. 

The  analysis  is  analogous  to  that  of  Doty  [12]  and  Doty  and  Rasmussen  [13] 
for  obtaining  approximations  for  hypersonic  flow  past  circular  cones  at 


small  angle  of  attack,  which  were  shown  to  be  very  accurate.  The 
analysis  is  cast  in  the  form  of  hypersonic  similarity  theory,  and  the 
results  are  presented  in  appropriate  similarity  form. 

The  perturbation  expansions  involve  a small  parameter  for  the  angle 
of  attack  and  for  the  measure  of  the  eccentricity  of  the  ellipse  cross 
section.  These  expansions  are  not  uniformly  valid  in  a thin  vortical 
layer  adjacent  to  the  cone  surface.  It  can  be  shown  [14,15],  however, 
that  at  least  the  pressure  and  azimuthal  velocity  component  are  valid 
across  the  vortical  layer.  Since  these  two  variables  are  of  most  importance, 
further  consideration  of  the  vortical  layer  will  not  be  undertaken. 


L 
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2.  UNYAWED  ELLIPTIC  CONE  BODY  AND  SHOCK  GEOMETRY 
2.1  Body  Geometry 

In  rectangular  Cartesian  coordinates,  as  shown  in  Fig.  2.1,  the 

unyawed  elliptic  cone  body  is  represented  by 

2 2 

-TT+^T=1  * (2-1} 

aV  b z* 

where  a = tan  0 and  b = tan  e.  are  the  tangents  of  the  semi  vertex 

d D 

angles  of  the  semiminor  and  semimajor  axes  of  the  elliptic  cone.  In 
terms  of  spherical  polar  coordinates,  also  shown  in  Fig.  2.1,  equation 
(2.1)  can  be  rewritten  as 

tan  em 


tan  e = 


where 


tan  em  5 


/ 1 + e cos  24 > 
_ /2  ab 


/a*  + bc 


= b/l-e  = a/l+e 


= /ib  (l-e2)1/4 


(2.2) 


(2.3a) 


e - b2  - a2 


(2.3b) 


The  parameter  e is  a measure  of  the  eccentricity  of  the  elliptic  cone. 

The  first  three  terms  of  the  Fourier-series  representation  of  equation 
(2.2)  are  given  by 

tan  0 = tan  ©m  [A  + cos  2<f>  + A^  cos  4<j>  + ...  3 , (2.4) 

where 


A„  --  - 


2 * *k2  /T+e 


A = 1 KU1 
0 1t  /T+e 

4 [(2-k2)  K(k)  - 2E(k)] 


(2.5a) 

(2.5b) 


Figure  2.1  Cone  Coordinates  and  Geometry 
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A = - — A - — A 
A4  ~ 3 o 3e  A2 


,/2  _ 2e 
K = ’PT? 


71/ 

E(k)=  j 


K(k)  = 


^1  - k2sin2u  du 


0 /l  - k2sin2u 


(2.5c) 


(2.5d) 


(2.5e) 


(2.5f) 


are  the  complete  elliptic  integrals  of  the  first  and  second  kinds.  For 
small  values  of  the  eccentricity,  e,  the  Fourier  coefficients  can  be  ex- 


panded as 

Ao  - 1 + 16  e2  + °<e4>  < 

A2  = ' I [ 1 + He2  + 0(e4^  ( 

A4  = fg  e2  [ 1 +§e2  + 0(e4)]  ( 

When  the  eccentricity,  e,  is  small  the  successive  Fourier  coefficients 


(2.6a) 

(2.6b) 


(2.6c) 


become  smaller  and  smaller.  For  sufficiently  small  e,  the  Fourier  coeffi- 
cients for  the  higher  harmonics  can  be  neglected. 

We  can  expand  equation  (2.4)  for  e in  a Fourier  expansion  and  ob- 


tain for  small  eccentricities 


where 


e = <5  - e COS  2$  + 0(e*') 

A A 2 

6 = tan_1(A0  tan  9m)  - sin  2 9m  + 0(e4) 

= 9m  + %K  [3  - 2 sin2em]  sin  29m  + 0(e4) 
m 32  L m m ' 


(2.7) 


(2.8a) 


A,  tan  9 A0£  0 A,4-  on  n 

' 5 ' H-A2  tan2e  CH^  * A«)  Si"  6"  + T 0+3flo  > si"  em  + >’ 
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e = | [1  + e2^  - | sin2em  + l sin4em>  + 0(e4)]  sin  2 0n 


(2.8b) 


The  parameter  e is  a new  measure  of  the  eccentricity  and  is  the  appropriate 
parameter  to  be  used  in  the  subsequent  analysis.  The  parameter  6 specifies 
the  semivertex  angle  of  the  basic  circular  cone  about  which  a perturbation 
analysis  is  to  be  performed. 

Comparison  of  the  two-term  approximation  (2.7)  with  the  exact  equa- 
tion (2.2)  for  the  elliptic  cone  is  shown  in  Fig.  2.2.  When  a = 0.2555 

and  b 3 0.3562,  such  that  e * 0.320  and  0 =16.30°,  the  two-term  approx- 

m 

imation  gives  a good  representation  of  the  actual  ellipse.  When  a = 0.2256 
and  b = 0.4034,  such  that  e = 0.523  and  = 15.56°,  the  representation 
is  not  as  good,  but  still  a reasonble  approximation  when  precise  accuracy 
is  not  paramount. 

2.2  Shock  Geometry 

The  conical  shock  wave  attached  to  the  elliptic  cone  is  assumed  to 


have  the  form  (for  a = 0) 


es  = B - eg  cos  2<j>  + 0(e  ) , 


(2.9) 


which  is  analogous  to  equation  (2.7)  for  the  body  shape.  Here  6 is  the 
semivertex  angle  of  the  basic  circular  shock  corresponding  to  the  basic 
body  with  semivertex  angle  S.  The  factor  g is  to  be  determined  from  the 
perturbation  analysis.  It_,in  effect, represents  the  deviation  of  the  shock 
eccentricity  from  the  body  eccentricity. 

By  means  of  vector  analysis,  we  find  the  unit  outward  normal  on  the 


shock  to  be  given  by 


?a4fn¥A'  % * °u2> 


(2.10) 


This  result  is  needed  to  establish  the  shock  jump  conditions. 
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3.  BOUNDARY  CONDITIONS 


3.1  Expansions  for  the  Flow  Variables 


Let  the  velocity  vector  for  conical  flow  be  represented  in  spherical 


coordinates  by 


V = u(e  »<t>)er  + v(  e ,4>  )eQ  + w(e,*)e  . 


(3.1) 


The  Fourier  representation  for  the  body  and  shock  shapes  suggest  that  the 
velocity  components,  pressure,  and  density  can  be  expanded  in  the  follow- 
ing forms,  valid  outside  the  vortical  layer  adjacent  to  the  body  surface 


(for  a = 0)' 


u(e,<j>)  * uQ(e)  + eu.j(9)  cos  2*  + 0(e  ; 

v(9,4>)  55  VQ(e)  + ev-j  (9)  cos  2$  + 0(e2) 

2 

w(9,4>)  = ew^ ( 9 ) sin  241  + 0(e  ) 

p(e,4»)  * PQ ( 9 ) + ep^e)  cos  24  + 0(e2) 

p(9,4>)  = Pq(0)  + epi  (0)  COS  2<j)  + 0(e2) 


(3.2a) 
(3.2b) 
(3.2c) 
(3. 2d) 
(3.2e) 


The  lowest-order  terms,  with  the  subscript  naught,  pertain  to  the  basic 
circular-cone  solution,  which  is  presumed  known. 


3.2  Free  Stream  Normal  Velocity  at  the  Shock 

The  free  stream  velocity  in  spherical  coordinates  is  given  by 


V = V cos  9 e - sin  9 e„ 

00  ao  y 0 


(3.3) 


At  the  shock,  9 = 9$  = e - eg  cos  241,  we  have  to  first  order  from 
cos  (eg  cos  2$)  £ 1 , sin  (eg  cos  24>)  £eg  cos  24> 


V = V {cos  e + eg  sin  3 cos  24>}  e 


* 2 1 

{sin  S - eg  cos  S cos  24>)  e + 0(e  ) (3.4) 

0 __ 
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Ustng  equation  (2.10),  we  find  the  normal  velocity  at  the  shock  to  be, 
to  first  order. 


V *n  = - V sin  b + egV  cos  b cos  2<t>  + 0(e  ) 

00  S 00  00 


(3.5) 


3.3  Downstream  Pressure  at  the  Shock 

The  shock  jump  conditions  give  for  the  pressure  ratio  across  the 
shock  the  well-known  relation 

Pc 


— = 1 + (M2  - 1) 

Poo  y+1  n 


(3.6) 


for  a thermally  and  calorically  perfect  gas.  Here  y is  the  ratio  of 


specific  heats,  y = c /c  , and  M = ^ *n/a  is  the  normal  free  stream 

P v n 


Mach  number,  a being  the  free  stream  speed  of  sound.  Substituting  (3.5) 

oo  v 

into  (3.6)  yields 

Pc 


= 1 + fvT  (Kg  - !)  -e  K2g  cot  6 cos  2 <P  + 0 (e2)  (3.7) 


where  K = M sin  B. 

B 00 


Equation  (3. 2d)  for  the  pressure  evaluated  at  the  shock  reads 

2, 


P(0S»4>)  = P0(9S)  + eP-|(0s)  cos  24>  + °(£t)  (3.8) 

Transferring  this  value  to  the  basic  unperturbed  shock  by  means  of  a 


Taylor  expansion  yields,  to  first  order, 

p(0$,<j>)  = PQ( 8)  + e |~-  (jjp  g + Pl(B)]  cos  2<i>  + 0(e2)  (3.9) 


Identifying  the  first-order  perturbation  terms  in  (3.7)  and  (3.9)  leads 
to  the  shock  boundary  conditions  for  the  perturbation  pressure 


^ . 4x  -2 


Y+,  9k;  cot  6 + i 


dPc 


de  J 


(3.10) 


B 


For  the  undisturbed  conical  flow,  the  pressure  gradient  can  be  evaluated 
in  terms  of  the  velocity  gradient.  Hence  we  have 


i 


? 


|—  j-dV  — r 

" 6 |_v. '°s  8 * IhtJ  J <3-’'> 


apoi  r 

-t-2-  = p V sin  6 V cos  8 + 

d0  j0  » - L - 


Thus,  alternatively,  the  pressure  boundary  condition  can  be  written 


P^S)  * p.9  sin  e coss  [j£y  + { 3r]  / v.cosB}  ] (3-12) 


The  derivative  (dv  /de ) can  be  evaluated  in  terms  of  the  density  ratio 

0 p 


across  the  shock,  as  shown  in  the  next  section. 


3.4  Downstream  Density  at  the  Shock 


The  pressure  ratio  across  the  shock  is  given  by 


p (y-i ) fT  + 2 

oo  n 


(3.13) 


(y+1) 


Expanding  this  expression  analogously  to  that  for  the  pressure,  we  obtain 


the  shock  density  to  first  order: 


P„  2 

— » £ + e£.  COS  2$  + 0(e  ) 

p_  o 


(3.14a) 


where 


’o  ' pQ(s) 


(y-i)kj  + 2 


(y+1)  k; 


(3.14b) 


5,  5 2g  cote[,0  - 


(3.14c) 


explicit  expression  for  p^(8)  Is  not  needed  It  can  be  shown  from  the 


basic  cone  solution  that 


£M-v«.[*-SU] 


(3.15) 


3.5  Downstream  Velocity  at  the  Shock 


The  velocity  components  immediately  downstream  of  the  shock  are  ob- 


tained from  the  governing  conditions  on  the  normal  and  tangential  compo- 


nents. The  normal  components  of  the  velocity  is  governed  by  the  mass  con- 


servation equation: 
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(V-n)$  = (V^-r^)  (3.16) 

The  rtght  side  is  determined  to  first  order  by  means  of  (3.5)  and 
(3.14): 

(V*n)  = - c V sin  B + eV  U g cos  B - £.  sin  B)  cos  2$  + 0(e  ) 

5 0 00  00  0 I 

(3.17) 


The  left  side  of  the  previous  two  equations  is  determined  to  first  order  from 
(2.10).  (3.1),  and  (3.2)  as 

(v*n)$  = VQ(es)  + £v1(es)  cos  2<t>  + 0(e2)  (3.18) 

Transferring  this  value  to  the  unperturbed  shock  yields 

( V - n ) $ = vo(b)  + e j^v1  (b)  - g ] cos  2$  + 0(e2)  (3.19) 

Comparing  equation  (3.17)  and  (3.19)  leads  to  the  following  results: 


vo(c)  = " V-  S1'n  3 


Vb)  = - C-jV^sin  B + g [50Vcos  B+  (^j  J 


Utilizing  (3.14c)  for  ^ now  gives 


v^B)  = g 


2(Y-1 ) 
Y+1 


’0  J 


V COS  B + 


rdv  -r 

b2)  ] 


(3.20) 

(3.21a) 


(3.2ib) 


The  other  two  velocity  components  at  the  shock  are  determined  from 

conservation  of  the  tangential  components: 

V x n = V x n (3.22) 

00  s s s 

Substitution  of  the  first-order  expansions  and  transfer  of  the  conditions 
to  the  basic  unperturbed  shock  yields 


u„(b)  = V cos  b 
0 00 

u,(b)  = g V sin  B(l-cJ 

I 00  0 


(3.23) 

(3.24) 
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w,  (b)  = - 2gV  0 ~0  (3.25) 

I 00  0 

The  needed  boundary  conditions  at  the  shock  have  now  been  specified. 

3.6  Surface  Boundary  Condition 

At  the  body  surface,  9 = ec  i 6 - ecos  2<j>,  the  normal  velocity  must 

vanish: 

( V - n ) c = 0 (3.26) 

By  means  of  vector  analysis,  we  find  the  unit  outward  normal  vector  on 
the  elliptic  cone  surface  to  be,  to  first  order, 


2 e sin  2<t> 


n = e " T ' T e + 0(e  ) 

c 9 sin  6 <j>  v ' 


(3.27) 


Substituting  this  expression  and  the  velocity  expansions  into  (3.25)  and 
transferring  the  boundary  conditions  to  the  basic  circular  cone  surface 
leads  to  the  surface  boundary  conditions 

(3.28) 


VQ( 6 ) = 0 


V-j  ( 5 ) 


dv  -i 
o 1 

de 


(3.29) 


Equation  (3.28),  of  course,  is  the  tangency  condition  for  the  basic  cone 
problem,  and  equation  (3.29)  is  the  surface  boundary  condition  for  the 
first-order  perturbation.  Rigorously,  we  should  have  obtained  (3.29)  by 
matching  the  outer  expansion  with  an  inner  expansion  for  the  vortical 
layer  adjacent  to  the  cone  surface.  Such  an  analysis  shows  that  equation 
(3.29)  is  indeed  proper  [14,15]. 
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4.  PERTURBATION  EQUATIONS 

The  pressure,  density,  and  velocity  are  governed  by  the  equations  of 
change  for  mass,  momentum,  and  energy,  plus  appropriate  equations  of  state. 
Here  we  assume  the  flow  is  inviscid,  nonconducting,  steady,  and  behaves 
as  a thermally  and  calorically  perfect  gas.  We  wish  first  to  obtain  ex- 
pressions for  the  pressure  and  density  perturbations  in  terms  of  the 
velocities,  and  then  finally  to  obtain  a single  equation  for  only  one  of 
the  perturbation  velocity  components. 

4.1  Energy  Equation 

For  steady  flow,  the  energy  equation  can  be  expressed  in  terms  of 
the  entropy,  s,  as 

^.grad  s = 0 (4.1 ) 

Thus  the  flow  on  either  side  of  the  shock  is  isentropic,  that  is,  the 
entropy  along  a streamline  is  a constant.  Of  course,  the  uniform  flow 
upstream  of  the  shock  is  homentropic  since  the  entropy  is  the  same  on 
every  streamline  before  it  passes  through  the  shock.  If  the  entropy  is 
expanded  in  the  form 

s(e,$)  = SQ(e)  + eS.|(e)  cos  2<p  + 0(e  ) (4.2) 

then  expansion  of  equation  (4.1)  leads  to  the  result  that  sQ  is  a constant 
and  that  s-j  is  also  a constant.  Thus  the  zeroth-order  flow  past  the  basic 
cone  is  homentropic  downstream  of  the  shock.  Since  s^  is  a constant,  the 
first-order  entropy  perturbation  depends  only  on  the  azimuthal  angle  !>. 

This  result  is  not  valid  on  the  body  surface  because  the  body  surface  is 
a stream  surface  that  has  a constant  entropy,  since  grad  s is  a vector  that 
is  perpendicular  to  the  surface  s ( 9 ,4> ) = constant  and  is  also  perpendicular 
to  the  velocity  V which  is  tangent  to  the  body  surface.  Thus  expansion 
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(4.2)  Is  not  valid  in  a thin  layer,  called  the  vortical  layer,  adjacent 
to  the  body  surface.  Only  the  pressure  and  azimuthal  velocity  are  uni- 
formly valid  across  the  vortical  layer. 

From  the  thermodynamic  state  relation  for  a perfect  gas 

(4.3) 


ds  - 4R  . Y dp 
S P p 


we  can  deduce  that  the  first-order  pressure  and  density  perturbation 
satisfy  the  relation 

~ M _ 

’1 


P-J  (©  ) 


P 1 ( 9 ) Sn 


P^eT  "Y^y 


= constant 


(4.4) 


O'  ' V 

Another  relation  between  the  pressure  and  density  perturbation  can 
be  found  from  energy  considerations.  Since  the  uniform  flow  upstream  of 
the  shock  is  homenergic  (constant  total  enthalpy)  and  homenergic  across 
the  shock,  it  is  also  homenergic  downstream  of  the  shock.  Thus  we  have 
for  a perfect  gas 

..  . „2  . p_  V2 

(4.5) 


„2  p r 

v D V V ~ c°  oo 

— 4-  ^ ^ — + -a-  = constant 

Y-l  p 2 y-1  2 


Substituting  the  perturbation  forms  (3.2)  into  this  equation  and  extract- 
ing the  first-order  perturbation  yields 

y-l  P' 


P1  (e)  p-j  (e) 


(u  U,  + V V . 

P0  0 1 0 1 


(4.6) 


This  is  a relation  between  the  pressure  and  density  perturbations  in  terms 
of  the  velocity  perturbation. 

Equations  (4.4)  and  (4.6)  can  be  solved  separately  for  the  pressure 
and  density  perturbations.  We  obtain 


P-j  ( 0 ) 


- p (u  U.  + V V ) 
0 0 1 0 1 


V l 


(4.7) 
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P^0)  = - PQ(uoUl  + vqVi)  + PqF1 


(4.8) 


where  aQ(0)  = yPq/po  is  the  speed  of  sound  squared  in  the  basic  unper- 
turbed flow,  and 


= constant 


(4.9) 


The  pressure  and  density  perturbations  are  thus  determined  in  terms  of 
the  velocity  perturbations.  The  constant  can  be  obtained  from  the 
shock  boundary  conditions.  Thus  we  obtain 

P-i  ( 6)  p (0)  — -r 

Fi  = L u°(6)  ul<6)  *,o(s) 


yg  V2  sin  S cos  6 (1  -S  )2/a2(s) 
°°  0 0 


(4.10) 


4.2  Momentum  Equation 

The  momentum  equation  for  fnviscid  steady  flow  is 


nr 


- y x u = - vp  , 


(4.11) 


where  n = curl  V is  the  vorticity  vector.  For  the  vector  quantities  let 


V = VQ  + evi  + 0(e) 


S2  = efi|  + 0(e  ) 


(4.12a) 


(4.12b) 


where  = curl  v^  and  the  spherical  components  of  the  first-order  pertur- 
bation velocity  v^  can  be  ascertained  from  equations  (3.2).  When  the 
perturbation  expansions  are  substituted  into  equation  (4.11),  the  first- 


order  perturbation  equation  can  be  extracted,  and  it  is 


i 

L 
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(V  \ 

(p^COS  2 <t>)  V + Pq^C^o*^)  - P0^0  X Q1  = - v ( p 1 cos  2$) 

(4.13) 

The  first-order  vorticity  perturbation  has  the  components 


“i  -ritr?  [a?  <V1n  8>  + 2vi  ] (4-’4a> 


= . sinja_  r2u  + w 

r sin  e [_/ul 


sin  e 


(4.14b) 


We  also  have 


cos2j> 


r ^ 

_ V1  " d0 


1 

J 


(4.14c) 


Vo  x n1  = er(Vl  ] - e0(uonl  ] + %(Vl  ’ Vl  } (4J5) 

* 0 r 


From  equation  (4.13)  it  can  be  determined  that 


(*o  X V*er  = 0 


(4.16) 


It  follows  from  (4.15)  that  n..  =0,  and  hence  that 

4> 

du. 


'I  ' de 


Thus  equation  (4.15)  reduces  to 


V„  x n,  = e (u  n-  - v fi.  ) 
o 1 4>  o 10  o lr 


The  ^-component  of  equation  (4.13)  can  be  written  as 


(4.17) 


(4.18) 


Vi„  - Vir  --  - 2 [ £ 4 v,  + Vl  ] rM4?  <4-,9> 


The  pressure  perturbation  p^  can  be  eliminated  by  means  of  equation  (4.8), 
and  we  get 


Vi9  ‘ Vir 


- 2a?  F 

o hl  sin  2<t 


r sin  6 


(4.20) 


where  a^(e)  = yPq/p0  and  F-|  is  a constant  given  by  (4.10).  By  means  of 
(4.14a,b)  we  can  further  rewrite  (4.20)  as 


2F 

*o  3?  [2ul  + "ls,n  9J  * uo  [2ul  + "lsi"  e]  * V-  ao(s) 


(4.21) 


Introducing  the  integration  factor 


9 

1(9)  s exp  r ( f de  ~[  , (4.22) 


we  can  integrate  equation  (4.21)  and  obtain 


2F  , / a^  I 

2u.  + w.sin  0 = — - r — — de 

11  y 1 J vo 

8 


(4.23) 


The  constant  of  integration  vanishes  by  virtue  of  the  shock  boundary  condi- 
tions (3.24)  and  (3.25). 

Equations  (4.17)  and  (4.23)  give  v1  and  w1  in  terms  of  u1 . It  remains 
to  find  a single  equation  for  u.j,' 

4.3  Continuity  Equation 


Mass  conservation  is  described  by  the  continuity  equation 

div(P\b  = 0 


(4.24) 


When  the  perturbation  expansions  are  substituted  into  this  equation  and 


the  first-order  perturbation  extracted,  we  get 


div  (p-|COs  2<f>  + PqV-j  ) = 0 


(4.25) 


With  the  basic-flow  result  div(pQVo)  = 0 accounted  for,  this  equation 
can  be  rearranged  to  read 


_ rp, cos  2<h  v. 

div  v,  = - V • v *Vp 

1 0 p0  p0  0 


(4.26) 
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If  it  were  assumed  that  the  density  varies  so  slowly  that  it  is  approx- 
imately a constant,  then  the  right  side  of  equation  (4.26)  can  be  set 
equal  to  zero  as  an  approximation,  and  a relatively  simple  equation  for 


v.j  ensues.  (The  right  side  does  vanish  at  the  body  surface,  0=  6 .) 


This  so-called  constant-density  approximation  leads  to  remarkably 
accurate  results  for  the  basic-flow  solution  and  for  the  flow  past  a 
circular  cone  at  angle  of  attack  [12,13].  These  constant-density  approx- 
imations were  justified  a posteriori  by  comparison  with  extensive  numer- 
ical tabulated  results.  Such  comparisons  are  not  so  available  for  the 
elliptic-cone  problem  of  present  interest,  and  hence  we  must  devote 
more  attention  to  the  terms  on  the  right  side  of  equation  (4.26). 

Tn  equation  (4.26),  let  us  now  eliminate  p-|(e)  by  equation  (4.7) 
and  utilize  the  following  basic-flow  results: 


2 2-^2 

a*"  dp  , dp  , da^  ,-V 

pq  de  y-1  de 


p0  d0 


d_ 

de 


(4.27) 


Equation  (4.26)  can  now  be  expanded  and  written  in  the  form 
-1  dv 

- + cot  e ! 1 - b I v,  + 

l 


[’ 


1 


de 


sin  e 


U1  = " 


2(2u^+w1  sin  e) 


T 
sin  e 


where 


A(e)  = V° 


"7 

ao 


(4.28) 

(4.29a) 


B(e)  = tan  e — j- 
ao 


dvo^ 

% + ar 


C(e)  = j 1 + (y-1)  — ^ 

a <■  a 


~z 

ao 


u.  + 


{ 2 + (y-1 ) | 

1 ao  J 

} 


dv 


de 


(4.29b) 

(4.29c) 
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The  factors  A,  B,  and  C in  equation  (4,28)  are  variable  coefficients, 
and  they  stem  from  the  right  side  of  equation  (4.26).  At  the  cone  sur- 
face e = <$o>  the  factors  A,  B,  and  C all  vanish  since  vQ( <5  ) = 0.  At 
the  shock  surfaces3  6,  the  factors  A,  B,  and  C take  the  values 

25, 


A(B)  = (y+1  ) - 50(y-D 


Bf0)  = (y+i  ) - O^TJ 


(4.30a) 


(4.30b) 


1 


C - 2 cot  8 

C(6)  = - 2C+3 , (4.30c) 

(y+1)  - 5q(y-1) 

where  the  undisturbed  shock  density  ratio,  5Q,is  given  by  (3.14b).  In 
the  hypersonic  limit,  K.  = ® and  A(b)  takes  the  value  (y-1 ) /2y . Thus 

P 

for  y = 7/5,  A( 8)  = 1/7  for  K0  = For  hypersonic  flow  it  is  thus  safe 

P 

to  neglect  A(e)  compared  to  unity.  Since  A(9q)  = 0,  it  is  thus  reason- 
able to  neglect  A(e)  always  except  perhaps  very  near  sonic  conditions, 
for  which  A(8)  = 1 when  K.  = 1 . Likewise  C(e)  is  small  compared  to 

P 

p 

(2-4  esc  e)  in  the  coefficient  of  u-j  for  all  values  of  K^.  Thus  C(e) 
should  be  negligible  in  general  since  it  also  vanishes  at  9 = 6 . The 
factor  B(e)  varies  between  (y+1)/y  for  K = °o  and  2 for  KQ  = 1 , and  thus 

P P 

apparently  should  not  be  neglected  compared  to  unity.  Moreover,  the 

factor  (1-B)  changes  sign  between  the  shock  and  the  body.  On  the  other 

hand,  if  A,  B,  and  C are  all  neglected,  the  resulting  differential 

equation  still  behaves  properly  at  the  body  surface  since  A,  B,  and  C 

vanish  there  anyway.  In  addition,  the  resulting  second-order  differential  equation 

must  satisfy  the  two  boundary  conditions  for  u^s)  and  v^b)  at  the  shock. 
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r 


Thus,  even  though  A,  B,  and  C be  neglected,  the  resulting  solution 
should  behave  properly  at  both  sides  of  the  shock  layer,  and  presumably 
also  in  between.  As  mentioned  earlier,  neglecting  A,  B,  and  C yields 
good  results  for  the  circular  cone  at  angle  of  attack. 

Even  though  we  intend  eventually  to  obtain  an  approximate  solu- 
tion by  neglecting  A,  B,  and  C,  we  can  justify  the  results  still  further 
by  recasting  the  differential  equation  as  an  integral  equation.  Replacing 
v.j  by  (4.17)  and  w-j  by  (4.23),  we  can  rewrite  equation  (4.28)  as 

2 4F1  Hn(9> 

u"  + cot  0 u-j  + (2  - 4csc  6 )u,  = 2~  + H,  ( 0 ) (4.31 ) 

Y sin  0 


where  0 2 . 

<1  1 


Ho(0)  = \ \ vo  d0 

(4.32a) 

0 

H1  (e)  = Au!|'  + B cot  @ uj  + Cu^ 

(4.32b) 

o 

When  Hq  and  vanish,  a solution  to  equation  (4.31)  is  u^  = esc  0. 
This  suggests  the  substitution 
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where 


X(B)  = u1 (g)  sin  0 


= g sin  6 (1  - eQ) 


X*  (s)  = g sin2e  [(^L  - 3co 


V cos  6 + 


dv  >1  -r 

*41 


When  X(e)  Is  known,  the  velocity  components  are  determined  by 


11,(0)  = X(e)/sin  e 


(4.36a) 


(4.36b) 


(4.37a) 


v-j ( 9 ) = jx1  (e)  - 2X(e)  cotej/sin2  e (4.37b) 

2F 

w,  (e)  = jj-  2u,(0)  + HQ(9)J  / Sin  0 (4.37c) 


The  shock  eccentricity  factor  g is  determined  by  satisfying  the  surface 
boundary  condition  (3.29)  and  then  solving  for  g.  Equation  (4.35)  is  an 
integral  equation  because  the  unknown  function  X(9)  also  appears  in  the 
function  H,(0). 
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5.  APPROXIMATE  SOLUTION  FOR 
HYPERSONIC  FLOW 


5.1  Basic  Cone  Flow  Approximations 

For  hypersonic  flow  in  the  limit  M^+^and  sins  -*-0  such  that  the 
combination  K = M sin  9 remains  finite,  the  basic  cone  flow  can  be  approx- 

oo 

imated  accurately  by  [10,11] 


«Ue) 


= l - 


. 2 

sin  6 


f slsil  ♦ fslni-S]  1 (5.1a) 

L-  sin  6 '■sin  dJ  -1 


v (0) 


= - si 


L sin  9 -1 


and  the  relation  between  the  shock  and  body  angles  is  given  by 

sin  B _ A 
sin  6 J 


y+T 


1 


(5.1b) 


(5.1c) 


(M^sin  6)‘ 

Let  us  now  restrict  ourselves  to  small  angles  such  that  sin  9 = 0 and 
neglect  second-order  terms  in  (5.1).  Then  we  get 

u 


1 


(5.2a) 


V ‘ 9 


52 

~z 

0 


(5.2b) 


6 / 2 1^2 


(5.2c) 


where  K,  = H i is  the  hypersonic  similarity  parameter. 

0 00 


5.2  Evaluation  of  the  Integrals 

With  the  approximations  (5.2a,b),  the  integrating  factor  I becomes 


i = 
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p2  - *2 

L7T7 


] 


1/2 


(5.3) 


In  the  evaluation  of  the  function  HQ(e)  given  by  (4.32a),  the  speed  of 
2 

sound  squared  aQ(e)  varies  only  by  a few  percent  across  the  shock  layer. 

2 

Hence  we  replace  it  by  aQ(0),  and  the  function  HQ(e)  becomes 


Ho(9)  - { 1 - Jt~7  1 


We  then  obtain 

0 


a?(s) 


2 v 


!7~7 

1 82  - 62 


- n + 


«/  B2-62 


(5.4) 


cos"1 - cos"1 

6 


1 Ff 

B 6 


M9)  -[  a£(s)  r A2  _ x2 

-2-y-  de  j de  - 0 


0 


4V„ 


8"  - 6 
1 T 


0-6 


5fl2  . 
6 9 


(5.5) 


25  /s2-52 


f -1  5 

-1  6 

«2  + 

52  2 1 

[cos  J 

"C0S  8j 

4.  4- 

6 

T " 6 J 

(5.6) 


5.3  Approximation  for  the  Velocity 

The  first  approximation  solution  for  the  integral  equation  (4.35) 
is  obtained  by  neglecting  the  integral  involving  H^(e),  which  is  the  same 
as  omitting  the  factors  A,  B,  and  C in  the  original  differential  equation. 
We  thus  get 
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25  /"?  ^ 


6 9" 


cos"1  - COS 


l TL 

7 

6 


B - 6 


■1  5]  62  62  ~[ 

9J  6?  303  J 


(5.8) 


f 1 62 


„i(9)  (,4.P4)  .f, 


26 /f2  ^ 


.-1  £ 
6 


6-5 
~2  52 


B - 6 


Hie  .6 

cos  - - cos  T-  + — 3 + — y 
9 6eJ  3eJ  -1 


(5.9) 


where 


- 2F.ad(a)  2 

f1  s L2 — = - sin  B cos  6(1  - 5Q) 

' ¥oo 

- - 29V„6(1  - if 


X(B)  * gv  jr(l  - cn) 


X'(6)  > 962  [ ($  - 3C0)  v.  * & ] 


2 2 
" W 9V»S 


(y-1)  - C, 


.(y+1)  1 


J 


(5.10) 

(5.11a) 

(5.11b) 

(5.11c) 


Expression  (5.11c)  is  the  approximation  consistent  with  the  velocity 
approximation  (5.26b),  that  is, 


fdv  -i 
_! 2.1 
jde  I 


6 


- V i + ± 

00 

(5.12a) 

e 

- V2  - H 

(5.12b) 
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slnce 


£ = 1 - V 

0 B2 


(5.12c) 


5.4  Shock  Eccentricity  Factor 


The  shock  eccentricity  factor,  g,  can  now  be  determined  from  the 


condition 


= ldrj5=-2V* 


(5.13) 


Solving  for  g from  equation  (5.8)  yields 


I = (1  - £ )2 

g o' 


Bcos"1  r 


2/b4  - 


3 ^ 3 

^3-3?  + 0 -£fl)Sg. 

663  35  > 0 63 


+ l/e  -nil  fi+£l 

2 \Co  y+1  i U 53 
Denoting  the  basic-cone  shock-body  ratio  by 


(5.14) 


we  then  get 


i-  i r J cos  l< 

9 ^ L /ttt 


- B _ /x+1  + 1 

a * 5 ‘ / 2 .2  * 


3 cos"1  [-)  c c 9 
a . 6 / 6 . Z 


(5.15) 


, 6 / 6 2\  ^4  2 r 

7+T  (a  + CT  ] + 3a  ‘ CT  - 5_ 


(5.16) 


The  shock  eccentricity  factor,  g,  is  plotted  in  Fig.  5.1  as  a 
function  of  for  y = 7/5.  For  -*■  0,  which  corresponds  to  the  limit 

of  linearized  theory,  the  eccentricity  factor  tends  to  zero,  g -*■  0,  that 
is,  the  shock  tends  to  a circular  Mach  cone.  For  the  limiting  hyper- 
sonic flow,  g approaches  the  asymptote  g = 0.955,  and  the  shock 

tends  to  embrace  the  elliptic  cone  body.  The  shock,  however,  is  always 

less  eccentric  than  the  body.  When  K.  -*•  * and  y 1 , then  » -*■  6 and 

0 


g -*■  1,  and  the  shock  embraces  the  body,  in  agreement  with  hypersonic 
Newtonian  theory.  The  angle-of-attack  shock  eccentricity  factor,  g, 
obtained  in  Section  7,  is  also  shown  in  Fig.  5.1  for  comparison. 


5.5  Shock-Layer  Velocity  Profiles 

The  first  approximation  for  the  velocity  u.| , given  by  equation  (5.7), 
is  plotted  in  Fig.  5.2.  The  hypersonic  similarity  form  gives  u^V^S  as 
a function  of  e/6,  K.,  and  y.  Because  the  thickness  of  the  shock  layer 
varies  as  a function  of  K^,  the  shock  layer  is  normalized  by  means  of  the 
variable 


The  body  surface  corresponds  to  6 = 0 and  the  shock  surface  to  9 = 1 . 

At  the  body  surface  u,/V  6 is  insensitive  to  variations  in  Kr,  having 

I 00  0 

approximately  the  value  unity.  At  the  shock  surface,  u-j/V^  is  quite 
sensitive  to  variations  in  K&.  In  the  hypersonic  limit  = °°,  u-j/V^S 
Increases  only  slightly  from  the  shock  to  the  body. 

The  velocity  perturbation  v./V^  is  shown  in  Fig.  5.3  as  a function 
of  e,  y = 1.4,  and  various  values  of  K..  The  variation  of  v,/V  across 
the  shock  layer  is  analogous  to  the  variation  of  u^/V^S,  except  that 
v-j/V^  is  equal  to  -2,  as  imposed  by  the  boundary  condition  (5.13). 

The  azimuthal  velocity  perturbation  w^V^  is  shown  in  Fig.  5.4  as 

*\j 

a function  of  9,  y = 1.4,  and  various  values  of  K_.  At  the  shock  sur- 

<5 

face,  w./V  increases  as  K.  increases.  For  K = 2,  the  variation  of 

i °°  0 0 

w^/V^  across  the  shock  layer  is  very  slight.  At  the  body  surface,  w-j/V^ 

decreases  as  increases.  This  is  shown  also  in  Fig.  5.5.  When 

-►  0,  w ■|(<5)/Voo  -*■  -2,  which  is  in  agreement  with  linearized  theory.  In 


KS 


I 


Figure  5.1  Shock  Eccentricity  Factors,  y = 1.4 


Figure  5.2.  Radial  Perturbation  Velocity ,y  = 1.4 
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the  hypersonic  limit  K.  ■+■  <«,  w , ( 6 )/ V becomes  asymptotic  to  the  value 
0.659  for  y = 1.4.  The  corresponding  angle-of-attack  contribution, 
8ft(«)/V.,  discussed  in  Section  7,  is  shown  in  Fig.  5.5  for  comparison. 

5.6  Evaluation  of  the  Approximate  Analysis 

Let  the  first  approximation  for  the  integral  equation  (4.35)  be 
denoted  by  uj°\  which  is  given  by  equation  (5.7).  Then  the  integral 
equation  (4.35)  can  be  written,  for  small  angles,  as 

= uj0)+  ujc)  , (5.18) 


where 


U(C)(6)  = — 

sin  e 


• 3 n 

sin  9 


H,(0)  > 

7fiTe  ds]  d9  <5-'9a> 


f9  3 f6  ) 

I 0 del 

b b 


H^(e)  = Au!j  + B cot  e u.j  + cu-j 


(3.19b) 


(5.20) 


Consider  now  an  approximate  evaluation  of  the  correction  function 
(cl 

u^  ' . Towards  this  end  we  utilize  the  well-known  convergent  iteration 

procedure  for  Voltera-type  integral  equations.  The  first  step  is  to 

substitute  the  first  approximation  uj1^  into  H^(e)  and  then  evaluate 

ujc).  Expression  (5.7)  for  uj0^  is  fairly  complicated  and  does  not 

(c) 

lend  itself  to  a simple  analytical  evaluation  of  u|  . Figure  5.2  shows, 
however,  that  uj°^(e)  can  be  approximated  by  the  simple  formula 

u(0)  , ,2 

rr-lf)  g<6'k5>  (5-2,a> 


where  G(e,K.)  is  a slowly  varying  function  of  e and  K.,  being  approxi- 

0 0 

mately  unity.  It  now  follows  approximately  that 

'(0)  „(0) 


" - 2 ? 6 


(5.21b) 
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UV(0>  53 

T — ■ 6 7 G 


(5.21c) 


since  G is  nearly  a constant.  Again  it  may  be  verified  that  expression 
(5.21b)  gives  a reasonable  approximation  to  Fig.  5.3. 


Substitution  of  expression  (5.21)  into  (5.20)  and  than  (5.19b) 


yields 


U<c)(e)  - ,{«(.)  [e3  [ f f S&.a+VUU 

1 ‘ > . L 9°  9 eJ-l  j 


(5.22) 


In  the  evaluation  of  A,  B,  and  C,  given  by  equation  (4.29),  we  set 
2 2 2 

aQ(e)  = a0(<$)  since  aQ(9)  varies  by  only  a few  percent  over  the  shock 
layer.  With  the  approximations  (5.2a,b)  we  now  get,  for  small  angles. 


and  with  z h 9/6, 


2 1 
A = Nz  1 - 


(5.23a) 


B = N (l  - ly)  [ 2 + (y-1  ) NZ2  fl  - (.23b) 


where 


2 z2  1 l2 

C = - (y-1)  N2  K 1 - '-j 

0 Z > 


a^B)  (c  -l)(2a  +Y-1) 


(5.23c) 


(5.24) 


where  a = B/6,  given  by  (5.15).  With  the  formulas  (5.23),  the  integral 
in  (5.22)  can  be  evaluated.  We  obtain 


first  approximation  is  about  one  percent  or  less.  As  K.  decreases,  the 

6 

correction  becomes  slightly  more  negative,  being  approximately  -0.015 

at  K,  = 1.7.  As  K.  decreases  further,  the  correction  factor  increases 

0 0 

from  the  minimum  negative  value  and  becomes  zero  at  = 1 . Further 
decreases  in  K.  give  a rapid  increase  in  the  correction  factor,  becoming 

<5 

0.75  in  the  limit  K6  = 0 which  corresponds  to  linearized  theory. 

In  the  hypersonic  flow  range,  we  can  expect  the  first  approximation, 
(c) 

obtained  by  neglecting  u^  , to  be  very  accurate.  In  fact,  this  is  true 


r • 
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when  Kp  £ 1.  When  K.  < 1 , the  first  approximation  is  less  accurate. 

6 6 

We  note,  however,  that  the  correct  limiting  results  for  K = 0,  which 
correspond  to  the  limiting  case  of  linearized  theory,  are  recovered. 
This  is  true  because  the  perturbation  u-j  itself  vanishes  when  S ■+  0 
and  also  because  the  surface  boundary  condition  v,(6)  = -2V  is  en- 
forced  on  the  first  approximation  itself  (which  produces  the  first 
approximation  for  g).  The  surface  boundary  condition  is  always  ex- 
actly enforced.  Thus,  although  the  linearized  limit  = 0 is  re- 
covered, the  approach  to  the  linearized  limit  is  in  error.  Because 
the  results  are  correct  at  K = 0 and  very  nearly  correct  at  K = 1 , 

the  error  in  the  range  0 < K.  < 1 , while  greater  than  for  K > 1 , 

6 <5 

is  less  than  that  indicated  in  Fig.  5.6. 

The  above  observations  are  born  out  by  the  results  of  Doty  and 
Rasmussen  [13]  for  hypersonic  flow  past  a circular  cone  at  angle  of 
attack.  Their  approximate  analysis  was  analogous  to  the  present  analy- 
sis, but  the  results  could  be  compared  extensively  with  well-known 
tabulated  results.  The  agreement  was  very  good.  Thus,  whereas  ex- 
tensive tabulated  results  for  the  elliptic  cone  do  not  exist  as  they 
do  for  the  circular  cone  at  angle  of  attack,  the  above  error  analysis 


and  the  experience  of  Doty  and  Rasmussen  justify  confidence  in  the 
present  analysis. 
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6.  SURFACE  PRESSURE  ON  THE  ELLIPTIC  CONE 

6.1  Surface  Perturbation  Pressure  Coefficient 

The  pressure  on  the  surface  of  the  cone  is  given  by  (for  a = 0) 

p( ec »‘P)  = PO(0C)  + eM0c^  cos  + (G-D 


where 


0c  = 6 - eCOS  2()>  + 0(e  ) 


(6.2) 


Transferring  the  expression  (6.1)  to  the  basic  cone  gives 

p( ©c  »<i> ) = P0U)  + e Tp ^ ( 6 ) - J COS  2$  + 0(c2) 

5 (6.3) 

At  the  cone  surface,  the  gradient  of  the  basic  pressure,  pQ,  vanishes 

since  v (<s)  = 0,  that  is, 
o 


dp. 


Id0  J 


- »0<5>  »0<5>  L“o  *3r  j ■ 0 <6'4> 


Hence  we  have 


P(ec>4>)  = P0(«)  + ep-j  (5)  COS  2<t>  + 0(e)  (6.5) 


The  pressure  coefficient,  Cp,  is  defined  by 

E--  1 

c - P~ 

P J M2 
2 « 


Thus  we  can  write 


where 


Cd  = Cd  + eCD  cos  2*  + 0(e) 
P P0  P] 


P«(fi) 


= P„ 


- 1 


0 i 


(6.6) 


(6.7) 


(6.8a) 


2p1C«) 

*P1  ’ ypM*’ 


(6.8b) 


The  pressure  perturbation  is  given  by  equation  (4.8),  and  since  vq( 6 ) = 0, 


we  have 


P-|(<5)  = - P0(«)  u0(6)  ^(5)  + P0(a)  F1 


(6.9) 


Using  the  previous  results  for  the  first  approximation,  we  obtain  the 


hypersonic  similarity  form 


'P1  2N  Po(5)  ra 
5 " 7 P.  L „3 


a _ a0(e) 


a‘(6) 


I 

o-  (6-10) 
00  — 


where 


p0<4)  - , + 1 k2  fio 

p 2 6 2 

roo  ^ 


(6.11) 


and  N is  given  by  (5.24).  From  an  analysis  of  the  basic  cone  flow 
[11 ,12,13]  we  have 


**> 


Vb) 


= 1 + 


/ 1 \ 2 i 21  -i 

lY-1 } a In  a + - 1 

a — 

(a2  - 1 )(2a2  + y - 1 ) 


(6.12a) 


P 2.2 

Ko  _ n , a In  a 

“3 1 + — 7^ 

6^  a - 1 


(6.12b) 


Figure  6.1  shows  C /$  plotted  as  a function  of  K.  for  y = 1.4. 

P1  r 5 

For  K.  = 0,  we  get  C /6  = -2,  in  agreement  with  linearized  theory.  In 

6 p-j 

the  hypersonic  limit,  K = °°,  C /6  approaches  the  value  -3.811.  Thus 

6 P] 

in  the  hypersonic  limit,  the  surface  pressure  coefficient  approaches 
the  form 


= 2.094  - 3.811  j cos  2*  + 0(e2) 


(6.13) 


I 


• ^ 


f 
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6.2  Comparison  with  Experiment 

Surface  pressures  were  measured  on  two  different  elliptic  cone 
models,  each  at  free  stream  Mach  numbers  3.09  and  6.0,  by  Zakkay  and 
Visich  Ci 6 1 - The  geometric  properties  of  these  models  were  as  follows: 


Model  I 
a = 0.2555 
b = 0.3562 
e = 0.320 


Model  II 
a = 0.2256 
b = 0.43®  .4P34- 
e = 0.523 


6 = 16.36°  0 = 15.56° 

m m 

e/6  = 0.155  e/6  = 0.266 

6 = 16.64°  =0.29 0+  rad  6 = 16.28°  = 0.2841  rad 

These  two  models  have  the  same  cross  sectional  areas  for  the  same  sta- 
tion along  the  elliptic-cone  axis.  A circular  cone  with  the  same  cross 
sectional  area  has  a semi-vertex  angle  of  16.79°.  The  experimental 
data  are  compared  with  the  results  of  the  present  analysis  and  also 
with  the  analysis  of  Martellucci  [5].  Martellucci  used  an  extension  of 
the  so-called  linearized  characteristics  method.  This  method  essen- 
tially uses  the  first-order  perturbation  equations  utilized  in  the 
present  analysis,  but  takes  a finite  number  of  terms  in  a Fourier 
series  expansion  to  represent  the  shape  of  this  surface.  Martellucci 
used  the  first  six  terms,  (In  view  of  the  fact  that  the  Fourier 
coefficients  for  an  ellipse  decrease  in  magnitude  like  powers  of  the 
eccentricity,  e,  as  seen  in  equation  (2.6),  it  would  seem  that  a higher- 
order  theory  should  be  utilized  to  accommodate  the  higher-order  Fourier 
coefficients.  In  this  sense,  the  lineari zed-characteristics  method 


V 
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does  not  seem  to  be  entirely  rational,  at  least  for  the  ellipse.)  The 
perturbation  equations  were  integrated  numerically  by  Martellucci, 
but  the  general  results  corresponding  to  the  present  analysis  were  not 
obtained. 

Figure  6.2  shows  the  pressure  distribution  on  one  quadrant  of 
model  I for  the  elliptic  cone  for  M^  = 6.0,  which  corresponds  to 
= 1.747.  The  present  results  agree  very  well  with  the  data  on  the 
major  and  minor  rays  ($  = 90°  and  180°),  but  give  pressures  that  are 
too  large  in  between.  The  overall  agreement  is  good  considering  that 
the  perturbation  theory  should  be  strictly  valid  when  s « 6,  and  this 
condition  is  met  only  marginally.  The  results  of  Martellucci  give 
slightly  bettei  agreement  with  the  data,  but  probably  not  enough  to 
justify  the  numerical  computational  effort. 

Figure  6.3  shows  pressure-distribution  data  on  model  II  for 
M^  = 6.0,  which  corresponds  to  = 1.724.  This  model  is  substan- 
tially more  eccentric  than  model  I,  and  the  condition  e « 6 is  cer- 
tainly not  satisfied.  Nevertheless,  fairly  good  agreement  with  the 
present  analysis  definitely  represents  the  overall  trends  of  the  data. 

The  results  of  Martellucci  do  not  appear  to  give  substantially  better 
overall  agreement  with  the  data. 

Figure  6.4  shows  the  pressure-distribution  data  on  model  I for 

( 1 

the  smaller  Mach  number,  M =3.09,  for  which  K.  = 0.900.  Again  the 

OO  Q 

present  results  give  good  agreement  with  the  data  on  the  semi -major 
ray,  but  the  data  are  lower  otherwise.  The  overall  agreement  does  not 
seem  to  be  quite  as  good  as  for  the  higher  Mach  number  M^  = 6.0, 

which  might  be  partially  expected  on  the  grounds  that  the  approximate 
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analysis  is  less  accurate  at  K.  = 0.900  than  for  K.  = 1.747.  The 

0 6 

results  of  Martellucci  give  a little  better  agreement  with  the  data 
between  the  major  and  minor  rays. 

Figure  6.5  shows  the  pressure-distribution  data  on  model  II  for 
M = 3.09,  for  which  = 0.888.  The  agreement  with  the  present  analy- 
ses is  fairly  good  near  the  major  and  minor  rays,  but  poor  in  between. 
Again,  for  this  large  value  of  eccentricity,  higher-order  perturbation 
terms  are  probably  required.  The  1 inearized-characteri sties  method 
used  by  Martellucci  picks  up  higher-order  harmonics  in  Fourier  repre- 
sentation of  the  ellipse  and  thus  shows  somewhat  better  agreement  with 
the  data  between  the  major  and  minor  rays.  In  view  of  the  large  value 
of  e,  however,  there  are  probably  higher-order  perturbation  terms  that 
are  of  the  same  order  of  magnitude  as  the  higher-order  Fourier  har- 
monics. 


6.3  Drag  on  the  Elliptic  Cone 

The  pressure  force  on  a finite-length  elliptic  cone  is  given  by 


f = - 


P ( 0 c ) n dS 


(6.14) 


where  S is  the  area  of  the  slant  face  , and 


9C  = 6 ■ e cos  2(j>  + 0(e  ) 


(6.15a) 


n 


e9  ' 


2e  sin  2* 
sin  6 


0(s2) 


dS  = r sin  e 

c 

tan  9 

= z — 

cos  ec 


dr  d<j>  + 0(e2) 
dz  d<j>  + 0(e2) 


(6.15b) 


(6.15c) 
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A Experimental  Results  of  Zakkay  and  Visich 
O Numerical  Results  of  Mortellucci 


<p  DEGREES 


Figure  6.2,  Pressure  Coefficient  on  Body  Surface 
Model  I,  M^  = 6,  y = 1 .4 


Figure  6.3.  Pressure  Coefficient  on  Body  Surface 
Model  II,  M^  = 6,  y = 1 .4 
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£>,  Experimental  Results  of  Zakkay  and  Visich 
O Numerical  Results  of  Mortellucci 


! 


When  the  integration  is  carried  out  in  the  range  0 £ $ £ 2w  and 
0 < z < H,  where  H is  the  length  of  the  cone  axis,  we  get  for  the 
cone  drag 

D h f-ez  = poU)A  + 0(e2)  , (6.16) 

where 

A = ttH2  tan2  5 (6.17) 

is  the  base  area  of  the  basic  cone  of  semi-vertex  angle  <5.  Thus  the 

2 

drag  on  the  elliptic  cone,  ignoring  terms  of  order  e , is  the  same  as 
the  drag  on  the  basic  cone  of  semi -vertex  angle  6 having  the  same  length. 


‘ 


7.  ELLIPTIC  CONE  AT  ANGLE  OF  ATTACK 


7.1  Superposition  of  Results 

Let  the  freestream  wind  be  inclined  to  the  axis  of  the  elliptic 
cone  such  that  the  inclination  angle,  a,  is  measured  in  the  x-z  plane 
shown  in  Fig.  2.1.  Let  the  coordinate  system  remain  fixed  to  the  body 
so  that  the  z-axis  is  the  axis  of  the  cone.  The  surface  boundary  condi- 
tions remain  unchanged.  The  shock  shape  can  now  be  represented  by  [12] 

'x*  2 2 

0$  = e + ag  cos  <j>  - eg  cos  2<j>  + 0(a  ,ae,c)  (7.1) 

where  g is  the  eccentricity  factor  associated  with  angle  of  attack.  We 
assume  a and  e to  be  small  and  of  the  same  order  of  magnitude.  Since 
the  first-order  perturbation  equations  and  boundary  conditions  are  linear, 
the  angle  of  attack  problem  for  the  circular  cone  can  be  superposed  with 
the  elliptic-cone  solution  at  zero  angle  of  attack.  The  subsequent  re- 
sults for  the  circular-cone  angle  of  attack  problem  were  obtained  by 
Doty  [12]  and  Doty  and  Rasmussen  [13]  by  an  analysis  very  similar  to  the 
foregoing  analysis  for  the  elliptic  cone.  The  angle-of-attack  shock 
eccentricity  factor  was  found  to  have  the  form 

3 ♦ r3 . wn  . r o(,2.1)v2-|  r„  * 


, 2,  I-,  , 4a  ~[  r / 2 _ .1/2  l"1  . r ^2  nvl/2~f 

5 - 2(l+a  ) M - | ct ( ct  -1)  In  a + (a  -1) 

(7.2) 

where  a is  given  by  (5.15)  as  before.  This  function  is  shown  in  Fig.  5.1. 
The  first-order  expansion  for  the  flow  variables  outside  the  vor- 


tical layer  can  be  written  as 
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SP,  2N  P0(B)  P„(J)  rn.J,  a|fs)  »,(«)-, 

0 0 

In  the  calculations  of  Doty  [12]  and  Doty  and  Rasmussen  [13],  the  ratios 
2 2 

Po($)/po(e)  and  aQ(B)/a  (fi)  were  set  equal  to  unity  consistent  with 
their  "constant-density"  approximation. 

The  angle  of  attack  perturbation  pressure  coefficient  given  by 
(7.9)  is  shown  in  Fig.  7.1  When  = 0,  the  limiting  result  of  linear- 
ized theory  is  recovered,  t / 6 = -4.  When  K = the  limiting  hyper- 

P]  6 

sonic  value  is  & /s  = -4.0836,  Near  K.  = 1,,  there  is  a small  dip  in 
P-]  6 

the  curve.  As  approaches  zero,  there  is  a small  overshoot  in  the 

curve  which  does  not  occur  in  the  exact  theory.  Over  the  range  of  K., 

□ 

the  value  of  & /6  does  not  differ  greatly  from  -4.  In  the  hypersonic 
P1 

limit,  =°°,we  can  write 


-4  = 2.094  - 4.084 

6 


cos  - 3.811 


cos  2<t> 


(7.10) 


correct  to  first  order  in  aande.  Expression  (7.10)  indicates  that  a/6 
and  e/5  should  be  sufficiently  small  in  order  for  the  perturbation  analysis 
to  be  valid. 


7.3  Comparison  with  Other  Results 

The  present  results  for  surface  pressure  at  = 6 and  angles  of 
attack  of  a = 5°  and  10°  are  shown  in  Figs.  7.2  and  7.3  for  models  I and  II. 

The  results  of  Martel! ucci  are  shown  for  comparison  together  with  the  experi- 
mental data  of  Zakkay  and  Visich  for  a = 10°,  which  is  a large  enough  angle 
of  attack  to  make  for  a demanding  comparison.  The  present  results  agree  well  with 
Martelluci  for  both  models  for  a = 5°,  the  worst  agreement  being  on  the  leeward  ray, 


-51- 

4>  = 0°.  For  a = 10°,  the  two  results  are  in  fair  agreement  for  the  small 
eccentricity  ellipse,  model  I in  Fig.  7.2,  except  near  the  leeward  ray 
of  the  elliptic  cone.  For  the  large-eccentricity  ellipse  at  a = 10°, 
model  II  in  Fig.  7.3,  agreement  between  the  two  results  is  good  only  near 
the  windward  ray.  For  model  II  at  a = 10°,  the  combined  large  values  of 
a/6  = 0.607  and  e/6  =0.270  render  the  first-order  perturbation  theory 
invalid,  especially  near  the  leeward  ray  where  the  separate  perturbations 
are  additive. 


Experimental  Results  of  Zakkay  and  Visich  for  a 
Numerical  Results  of  Martel lucci  for  a = 10° 
Numerical  Results  of  Martel lucci  for  a - 5° 
Present  Results  for  a = 10° 

Present  Results  for  a = 5® 


(j)  DEGREES 


Figure  7.3,  Comparison  of  Surface  Pressure  with  Martellucci's 
Results,  Model  II,  M = 6,  y = 1.4 
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7.4  Normal  Force  on  the  Elliptic  Cone 

The  normal -force  on  the  cone  is  found  to  be 

?*e  = Cw  a \ M2P  A + 0(a2,e2,ae)  (7.11) 

X N c °°  ® ' 

a 

* 2 2 

where  F is  defined  by  (6.14),  A s iH  tan  5 is  the  base  area  of  the  basic 
cone  of  semi -vertex  angle  S,  and 


'N 


CP 

2 6 


(7.12) 


is  the  a-  derivative  of  the  normal -force  coefficient. 
The  moment  about  the  cone  vertex  is 


M = - 


p(ec)  r x ncdS 


Evaluation  of  this  integral  gives 

" ' [v  2 M"P"HA]  * 0(“2,t2, 

where 
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3 6 


ea  J 


c = — c 

3 N 

a a 


% 

cp 
i pi 


(7.13) 


(7.14) 


(7.15) 


is  the  a-  derivative  of  the  moment  coefficient. 

For  completeness,  we  can  rewrite  the  drag  on  the  elliptic  cone, 
(6.16),  in  terms  of  an  axial -force  coefficient.  If  the  base  pressure 
is  reckoned  as  p , then  we  have 


\ MfP  A + (a2,e2,ae) 


where 


'A  2 


C.  = C 
A P, 
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"p  ' 

Ho 
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and  C is  given  by  (6.12b). 
po  r r 

P1  po 

The  ratios  -y—  and  -y  depend  only  on  Kg  and  y.  They  are  in- 

6 

sensitive  to  when  is  large.  Thus  when  is  large,  that  is,  for 

hypersonic  flow,  the  main  effect  of  the  cone  cross-section  shape 


is  determined  from  6.  For  slender  cones  (small  9 ) we  have  from  (2.3a) 
and  (2.8a)  the  alternative  forms 


6 = /aF  [1  + 0(e2)] 

5 = b /I -e  [1  + 0(e^)] 


The  cross-section  area  of  an  ellipse  is  proportional  to  the  product  ab. 

Hence  when  the  cross-section  area  is  held  fixed,  the  force  coefficients 

2 

are  independent  of  the  eccentricity  e when  terms  of  order  e are 
neglected. 

On  the  other  hand,  when  b is  held  fixed,  then  6 varies  with  the 
eccentricity  to  the  first  power.  The  normal  force  and  moment  are 
independent  of  & when  is  large,  but  CA  decreases  with  increasing  e 
when  b is  held  fixed.  Thus  the  lift-drag  ratio  increases  when  b is 
held  fixed  and  e increases. 
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8.  CONCLUDING  REMARKS 

General  flow  field  results  for  the  hypersonic  flow  past  an  elliptic 
cone  have  been  obtained.  The  results  are  valid  for  large  Mach  numbers 
and  small  stream  deflections  such  that  the  hypersonic  similarity  para- 
meter, = Moq6,  is  fixed  in  the  limiting  process.  The  results  are 
more  accurate  for  large  K{  (K.  > 1),  but  the  proper  linearized  theory 
result  is  recovered  when  Kr  -*■  0.  The  ellipse  eccentricity  factor,  e,  in 
the  analysis  must  be  small  in  the  strict  sense  that  e/5  <<  1.  Comparison 
with  experimental  results  indicates  that  e/5  need  not  be  very  much  less  than 

unity,  but  merely  moderately  less  than  unity,  for  acceptable  engineering  re- 
sults. In  addition,  the  angle  of  attack  should  be  such  that  a/5  <<  1 
in  the  strict  perturbation  sense.  An  important  feature  of  the  analysis 
is  that  the  basic  circular  cone  angle, 5,  has  been  well  defined  in  terms 
of  the  geometric  properties  of  the  elliptic  cone. 

The  methodology  of  this  analysis  can  be  extended  to  other  cross- 
section  shapes.  Each  term  in  a Fourier  expansion  of  the  cross-section 
shape  can  be  handled  in  an  analogous  manner  and  accurate,  approximate 
analytic  results  obtained.  Strictly  speaking,  however,  the  cross-sec- 
tion shapes  should  deviate  only  slightly  from  a circular  cone  in  order 
for  the  perturbation  analysis  to  be  valid.  Moreover,  successive  Fourier 
coefficients  should  not  decrease  in  powers  of  the  basic  expansion  para- 
meter, for  then  corresponding  terms  of  higher-order  perturbations  be- 
come equally  important. 
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